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We study the shape and growth rate of necks between sintered spheres with dissolution- 
precipitation dynamics in the reaction-Umited regime. We determine the critical shape that 
separates those initial neck shapes that can sinter from those that necessarily dissolve, as well 
as the asymptotic evolving shape of sinters far from the critical shape. We compare our results 
with past results for the asymptotic neck shape in closely related but more complicated models 
of surface dynamics; in particular we confirm a scaling conjecture, originally due to Kuczinsky. 
Finally, we consider the relevance of this problem to the diagenesis of sedimentary rocks and other 
applications. 



I. INTRODUCTION 

Sintering is a surface-tension driven phenomenon in 
which soHd particles packed below their melting tem- 
perature are consolidated via the growth of necks and 
subsequent shrinkage of pores between the particles. 
Quantitative understanding of sintering is crucial to 
study evolution of porous morphologies in industry and 
in nature - key properties of ceramic and metallic 
powders, like their strength and degree of compact- 
ness, are c ontrolled b y the sintering proces s dQ^rmaij , 
199(: : lHerri ng. 1950; K ingerv and Berei 119551 iKuczinskv . 
1949|) . Packed snow flakes and ice spheres may be bonded 



together by necks well below their melting tempera- 
ture - a phenomenon that is important for understand- 
ing snow avalanches a nd glacier flows ("Colbcck", "1998'; 
Hobbs and MasonL [l964: ..Kingenl Il96(i .Kuroiwa. 1961; 
Maeono and Ebinumal Il983|) . Sedimentary and other 



granular rocks also undergo sintering processes that affect 
their porosity and pe rmeability (Hay and Evange ■198& 
|Jurewic~ and Watsonl Il985t (Visseit Il999tl . 

Sintering processes typically exhibit several stages, 
identified by the degree of porosity of the material and 
the nature of the dominant surface tension, as well as 
the kinetic route. In the early stage of sintering, when 
particles are barely touching each other, the driving force 
is the surface tension between the solid particles and the 
surrounding pore space - be it filled with vapor, solu- 
tion, or vacuum. The chemical potential /isur(x) associ- 
ated with surface tension at a point x on the surface is 
given by /Zsur(x) = aVmHipi), where Vm is a molecular 
volume in the solid phase, H = (ki + K2)/2 is the mean 
between the two principal curvatures K1.2, and a is the 
surface tension, here taken to be isotropic for simplic- 
ity. Typically, the contact zone between the particles is 
highly concave and has a negative H, whereas the parti- 
cles themselves are convex with positive H . This differ- 
ence in H, and hence in the chemical potential jUsur(x), 
gives rise to a net flow of solid molecules from the pe- 
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riphery to the touching zone and to the formation and 
growth of necks between the particles. 

At later stages of the process, necks become thicker and 
variations in H along the surface are diminished. Most 
of the pores are now concentrated inside grains and on 
grain boundaries, and the thermodynamic driving force 
for further compactification is not cr, but rather the grain- 
boundary energy. The dominant kinetic mechanisms in 
these stages are associated mainly with bulk transport: 
grain -boundary diffusion, plastic flow and latt ice diffu- 
sion llaermanLlliM ISwinkels and Ashbvl Il98ll) . 

Here we focus on early stage sintering where the mass 
transfer is controlled solely by dissolution-precipitation 
through a surrounding solution. Such kinetics is domi- 
nant when the solid phase is in coexistence with vapor or 
solution, and the temperature is not too close to the melt- 
ing temperature, at which thermal energy is high enough 
to activate fast surface diffusion processes. In addition, 
we assume that thermodynamics (energy) and kinetics 
(diffusion) associated with grain boundaries can be ne- 
glected, as is the case, e.g., with amorphous materials. 
We will limit ourselves to cases in which the chemical re- 
action rate for dissolution-precipitation between the solid 
and solution is much slower than typical rates of molec- 
ular transport through the surrounding media. Under 
such reaction-limited dynamics spatial variations of /isoi, 
the chemical potential in the solution, can be neglected. 
Moreover, we will consider particles in solution under 
'open' conditions, in which the total amount of solid ma- 
terial is not conserved, but the solute concentration, and 
hence /Xgoi, is fixed in time. These assumptions are unre- 
alistic for industrial sintering of powders, in which proper 
inclusion of molecular transport and other considerations 
can be crucial. 

The prototypical geome try for ear l y stag e sintering, 
introduced by Kuczinsky ijKuczinskvl Il949() . and stud- 
ied extensively since then, is of a cylindrical neck 
between two touching spherical particles. Assum- 
ing a characteristic shape for the evolving ne c k pro - 
file, it was shown llHerrind IT95(it |K uczms k3 I1949D . 
that for reaction-limited dynamics, dominated solely by 
dissolution-precipitation, the neck thickens in time as t^^^ 
at short times. If diffusion through the vapor is much 
slower than the reaction rate, it was argued that the 
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neck thickens as t^^^ l)Hobbs and Masoil Il964j) . Other 
kinetic routes give rise as well to various power laws. 
During the last 50 years, many experiments aimed to test 
these theories, and power law growth rates of necks were 
indeed found llHobbs and MasonL Il964l: iKingervl Il960t 
iKingerv and Ber For sintering dominated by 

surface diffusion, several researchers argued, based on nu- 
merical simulations and analytic arguments, that Kuczin- 
sky's growth rates are not correct (Nichols and MuUins 
1965, German and Lathrop 1978, Eggers 1998). How- 
ever, the multitude of microscopic mechanisms for sinter- 
ing, and the fact that in real physical processes several 
of them may be significant, makes the verification of a 
simple theory a challenging task for the experimentalist. 

In this paper we present a detailed study of sintering 
driven by our model dynamics in this geometry. We focus 
on two aspects of the surface dynamics: the critical sur- 
face, which is actually an unstable fixed point of the ap- 
propriate surface dynamics, and the asymptotic growth 
profile. Our analysis of the critical surface provides us 
with conditions on initial contact geometries that can 
become persistent sinters, and enables us to determine 
growth rates at very early stages of the sintering process. 
For the asymptotic profile, we verify the t^^^ growth rate 
for the neck, derive analytically the correct neck profile, 
and explain how it selects the appropriate growth rate. 

The paper is organized as follows: In section 2 we intro- 
duce our model system and derive an equation of motion 
for cylindrically symmetric surface evolution. In section 
3 we study the structure of the unstable fixed point of 
this equation (the critical surface), and perform linear 
stability analysis of its dynamics. In section 4 we discuss 
the evolution of surfaces far from the critical surface, and 
show how the asymptotic dynamics can be described as 
a self-consistent solution of the equation of motion. In 
section 5 we conclude and suggest future directions. An 
appendix contains some calculational details. 



II. EQUATION OF SURFACE EVOLUTION 

Let us consider a solid in coexistence with its solution 
in a surrounding liquid (or vapor), and denote the inter- 
face between the two phases as F = x(u,w). We assume 
a first order single-component chemical reaction between 
the solid and solution, characterized by Kf - the dissolu- 
tion rate of a flat solid interface to a 'fresh' (unsaturated) 
solution. We assume that the surface tension a between 
the two phases can be considered as isotropic (which will 
be true well above the roughening transition or for an 
amorphous solid). The normal velocity of the surface 
M„(x) into the fluid region is given by the difference be- 
tween dissolution and precipitation rates 



Un{x) ^ -Kf [1~ e TT-j 



(2.1) 



where the precipitation rate is controlled by the Boltz- 
mann factor associated with the difference A/i(x) = 



/isur(x) — /isoi(x) between the chemical potentials of solid 
and dissolved molecules on the two sides of the interface 
r. The chemical potentials /L(sur(x) and /isoi(x) are given 
by 



(x) = 2iymCrH{-x) + ^fiat 

c(x) 



Msur 

^sol(x) = fcTlog 



Csat 



Mflat 



(2.2) 
(2.3) 



where /iflat is the chemical potential of a flat surface in 
equilibrium with a saturated solution at concentration 
Csat, c(x) is the concentration near the surface point x, 
Vm is the molecular volume in the solid, and an ideal 
solution is assumed. Here and elsewhere we define all 
concentrations with respect to the concentration in the 
solid. Equation (|2.2() follows from the geometrical iden- 
tity H{x) — {1/2){SS/SV), relating the mean curvature 
and the local variation in surface area SS with respect 
to a volume change SV of the solid. For our system we 
assume that the surface tension energy is much smaller 
than kT, and that the solution is nearly saturated. Equa- 
tion (|2.1|l then reduces to: 



w„(x) 



-K 



f- 



ff(x) 



kT , c(x, t) 
log 



Csat 



(2.4) 



In order to fully describe the evolution of a solid-fluid in- 
terface r(i), equation H2.4() must be supplemented by an 
advection-diffusion equation governing the concentration 
in the solution, for which c(x, t) is the boundary value at 
the interface with the solid phase. 

In this study we are interested in kinetic conditions 
for which any inhomogeneity in the concentration c(x, t) 
equilibrates much faster than the typical time for chem- 
ical reactions to cause significant changes of the sur- 
face shape by dissolution-precipitation events. In other 
words, we assume: 



T'diss S> Tdiff 



Tdiss Tadv ; 



where the kinetic time scales are defined by 
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(2.5) 



(2.6) 



where D is the diffusion constant, V is typical velocity 
of the advecting flow, Lpar is a typical size of a parti- 
cle, and Lporc is a typical separation between particles. 
Such kinetic conditions are typically realized in sedimen- 
tary rocks that dissolve extremely slowly, such that their 
surface morphology evolves over geological time scales. 

Under this reaction-limited dynamics, we can assume a 
homogeneous distribution of the concentration and hence 
the chemical potential in the solution: c(x, i) — > Coo(i), 
fixed by conditions outside the local region of interest. 
The surface dynamics equation (|2.4|l then reduces to: 



M„(x) 



-K 



'W 



f- 



[i?(x) - H4 . 



[kT/(2j^^a)]\og(c^{t)/c,,t) 



(2.7) 
(2.8) 
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Equation (|2.7|l takes the form: 



FIG. 1 Two sintering spheres-the sinter preserves the cylin- 
drical symmetry around the C^-axis. The width of the neck 
divided by the diameter of the sphere defines the small pa- 
rameter £. 



For = equation (|2.7|) becomes the well-known AUen- 
Cahn equation, describing the decay to global equilib- 
rium of a binary sy stem under the influence of surface 
tension l)Bravlll994|) . 

In this paper we focus on the evolution under (|2.7|l of 
a sinter between spheres of radius R. Typical curvatures 
in the neck region are much higher (in absolute value) 
than the curvature of the spheres, and therefore we take 
Coo{t) = Csat exp{2vmO'/kTR), which guarantees that be- 
fore contact is made the two solid spheres are in equilib- 
rium with the solution. This condition replaces the more 
common condition of conservation of solid matter, which 
is usually encountered in applications in material science, 
and which can be implemented in the open, reaction- 
limited case by suitably varying (slowly) with time. 

Depending on the surface geometry in the contact re- 
gion, the dynamics (|2.7|) may lead to dissolution of the 
contact, or to sintering - a neck growing between the 
two spheres. Evolution of surfaces under equation H2.7|l 
preserves cylindrical symmetry with respect to the axis 
connecting the centers of the two spheres. 



Let us rewrite equation H2.7() , taking advantage of this 
cylindrical symmetry. Defining the ^-axis as the axis of 
this symmetry, with p(C, t) the radial position of the sur- 
face (see Fig.^l, the in-plane (longitudinal) curvature k; 
and the out-of-plane (azimuthal) curvature 



+ " pTTTT^ 

yielding a mean curvature H of 

1. ^ 1 A P" 
H ^ -{Ka + Hi) = — , „ :; -2 



(2.9) 



(2.10) 



Furthermore, the normal velocity u„ is given by: 

Un^il+p'Y'^^dp/dt. (2.11) 

Upon rcscaling time and length by: 
t 



t 



kTl{R^KfV„,a) ' 



r = p/R, z = C/i?, 



dr 
di 



2\/l 



(2.13) 



In the following sections we will study the surface dy- 
namics described by equation H2.13|l . 



III. FIXED POINTS 

In this section we will discuss equation (|2.13(l near its 
fixed points - surfaces for which dr/dt = at each point. 
From equations H2.7I2.1U() we see that such surfaces are 
cylindrically symmetric surfaces of constant mean curva- 
ture (CSCMC) whose mean curvature H{z) ~ 1. First 
we discuss the structure of these surfaces, then we use 
linear stability analysis to study the dynamics near such 
surfaces. 



A. Structure 

According to equation (|2.13|) . CSCMC are defined by 
solutions to the equation: 



2VT" 



1 



r'2 = -- 



1 + r' 



(3.1) 



This is a second order ordinary differential equation 
(ODE) , and therefore two boundary conditions (BC) are 
required for a unique solution. Since we are interested 
here in sintering of two identical spheres, we consider the 
boundary conditions to be: 



r{z = 0) 



r'(z = 0) = 0, 



(3.2) 



where z = is the middle point between the centers of 
the touching spheres, e is the ratio between the neck size 
and the radius R of the unsintered spheres (which deter- 
mines iJ* ) , and the latter condition reflects the symmetry 
between the two sides of the contact. 



(2.12) 



Equations (|3.1I3.2|I can be easily integrated numeri- 
cally. The profiles r(z) and their derivatives r'{z) for two 
values of e are shown in figure 2, together with the profile 
rsphiz) of a sphere of radius 1, whose center lies on the 
2-axis at z = 1. Obviously, r'^ph[z) diverges as z ^ 0. 

The appearance of a singular inner region near z = 
in which |r'(z) — r'gpi^{z)\ >> 0(e) indicates that a reg- 
ular perturbative solution to equation H3.1|l is not possi- 
ble; this equation therefore calls for singular perturbation 
analysis. 



1. Singular perturbation analysis 

Another form of equation H3.1|l is achieved by consid- 
ering z as a function of r (assuming that the function 
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FIG. 2 Profile of a cylindrically symmetric constant mean 
curvature surface, obtained by numerical integration of equa- 
tions 13.113.21 . for two values of e: e = 0.1 (dots) , and e — 0.03 
(stars): (a) The profile r{z). The solid curve is a spherical 
profile Vsphiz). (b) The derivative of the profile r'(z). 



r{z) is one-to-one): 



r 1 + z'^ 
and the transformation of the BC (|3.2|l is: 



z{r = e) — 0, z' (r 



oo. 



(3.3) 



(3.4) 



Equations H3. 313.41) are more convenient for analytic 
study and will be used in the sequel instead of equa- 
tions H3.1I3.2|I . 

In order to construct a series expansion in e for z{r) 
in the interval r g one needs to consider sepa- 

rately the neck and the periphery, identify the small non- 
dimensional parameters in each of them, and then match 
the two expansions. 

Let us start with the peripheral region far from the 
neck, which we call the 'outer' region: r ^ e. In this 



region we can expand the profile Zout (f) as a series 

Zoutir)^z^^J,{r) + ezi]^,ir) + Oie'), (3.5) 

where z^'^^ir) is the profile of a sphere of radius 1, cen- 
tered on the z-axis at (r = 0, z = 1): 

z^'Kutir) = I - . 



(3.6) 



We now consider the next order in e. Substituting in 
()3.3|l the expansion (|3.5|l and the function z^^j (|3.6(l . we 
obtain the following equation for z^l^^: 

X'^out) ~ V'outI 



r(l — r^) 
whose solution satisfies 

(1) V _ 



(3.7) 



i^out) 



(3.8) 



where a is a constant of integration to be ultimately de- 
termined by a matching condition. Thus, to first order 
in e, the derivative of the outer solution is: 



>(1 - r2)3/2 



(3.9) 



A necessary condition for equation H3.5(l to be a valid 
asymptotic expansion is that: |(ZoJjj)'| ^ K-Zou*)'!- This 



leads to the conditions: 



r > 



1 



(3.10) 



Condition H3.10|l indicates the existence of an 'equato- 
rial' region 1 — ^/e < r < 1, where the outer pertur- 
bative expansion 1)3.5(1 is not valid, in which the profile 
is described by a different function z^q^ir) that matches 
Zoutir). As we shall see below, this actually amounts to 
a displacement of the singularity at r ^ = 1 in equa- 
tions H3.6I3.8|I to r ^ < 1. Since our primary interest 
here is the neck structure, we postpone the analysis of 
Zequ to the end of this section. 

The divergence of z' required by condition l|3.4|l indi- 
cates that a regular expansion like (|3.5|l cannot satisf y 
equation subject to the BC (E3I) fHinchl, Il99l|) . 

This fact implies that a perturbative solution to equa- 
tions (|3. 313. 4(1 must be constructed from an expansion 
whose 0*'* order term is different from (|3.6|) . We call this 
expansion Zi„(z), and write, similarly to equation H3.5|l : 



zUr)^zf^{r)+ez\^\r)+0{^) 



(3.11) 



The function 2^,1 (r) is required to satisfy equation 1)3. 3|l 
under the BC H3.4|l . and therefore we expect it to describe 
the profile near the contact between the two spheres. 

The BC (|3.4|l indicates that a natural length scale at 
the neck region is e, and therefore we re-scale, 



(3.12) 
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Equation H3.3|l then becomes: 



Assuming 



z' 


z" 


— 4 
r 


1 + [z'Y 


I' 




z" 


f 


1 


1 + {~z'Y 



(3.13) 



(3.14) 



and the expansion H3.11|l . we see that z^^^ = zf^^ /e satis- 
fies the equation: 



(5(0))' 



(g("))^' 
l + (i(o))'^ 



= 0, 



(3.15) 



which is just the equation of a catenoid (CSCMC with 
zero mean curvature). The solution of equation (|3.15|l 
subject to the BC (|3.4|) satisfies: 



(z("))' 



(3.16) 



Now let us calculate the first order term z^^ of the inner 
expansion. Substituting in equation (|3.13(l the expansion 
(|3.11ll and the catenoid function we get the following 



equation for z^^-* 



>le: 



?2 



whose solution satisfies: 

(z«)' = 



(3.18) 



The necessary condition |(z(^^)'| <C Kz*^"^)'! implies f <C 
1/Ve, or 



r < 



(3.19) 



which also establishes the validity of equation (|3.14|) . 
To first order in e, the inner expansion is thus 



1 



f2 



1 



1 



0(e2 



(3.20) 



From equations H3.10|l and (|3.19() we see that the inner 
and outer expansions are valid in non-overlapping regions 
in the interval [e, 1]. In order to match the two expansions 
Zin and Zouti one has to look for a separate solution z„^id^ 
valid in a 'middle' region around r ^ ^/e. Since the inner 
and outer expansions are valid at r <C ^/e and r 3> ^/^, 
respectively, a natural scaling for the variable and the 
function Zmid in the middle region is 



(3.21) 



The function Zmid must match, to first order in e, both 
Zin and Zout- This condition, and equations H3.9I3.2(J[I . 
imply the following form for the derivative z : 



2{f)=af + P +/(f)+0(e^ 
where asymptotically /(f) must satisfy: 



(3.22) 



\f{r)\<^f 
|/(f)| < af 



at f < 1 
at f > 1, 



(3.23) 



in order to enable matching to Zin and Zout ■ Substituting 
the form (|3.22(l in equation H3.3|l and expanding to 0{e) 
we arrive at the solution for /(f): 



1 ..3 

— e r 

2 ^ 



:-3 

- r 



3). 



(3.24) 



From equations H3.23I3.24|I we see that the asymptotic 
expansion H3.22|l for z^id is valid in the region 



< f < 



e < r < 1. 



(3.25) 



From equations (I;^.9l;^.l()l;^.19l;l:^()l;^.:^^l;^.:j4l;^:^ we 

conclude that to first order in e, Zmid{f) and Zin{r) match 
in the region e ^ r <C Ve, whereas Zmid{r) and Zout{r) 
match in the region ^/e <c r <C 1. The last matching 
condition also sets a = 1. 

To complete our first order asymptotic perturbation 
analysis we have to find z^qu{r) for r G (1 — y^, 1]. 
We note that unlike the BC 1)3. 4|l . which gives rise to 
a new singularity on top of the spherical profile Zoutii") 
in the neck region, the divergence of the outer expansion 
(|3.5|) for 1 — r ~ e, is caused simply because the origi- 
nal singularity, z'{r) ^ oo at = 1, is shifted to some 
rs(e) < 1. Analysis of equation H3.3|l around the singular 
point shows that Zequ{f) satisfies: 



equ 



c(e)/\/'"s(e) - r, 



1 - r < 1, 



where 



c(e) = l + 0(e), r,(e) = l + 0(e). 



(3.26) 



(3.27) 



The functions c(e) and 7's(e) are determined from match- 
ing Zequ to Zout to Oit) in the interval \/e <C 1 — r ^ 1. 
Performing the matching procedure yields 



1 



"equ 



r) 



(3.28) 



To conclude, to first order in e, the derivative z'{r) in 
the interval (e, 1 — e) is given by four different functional 
forms, equations H3. 913. 201 13.22l3.25|l , each valid in a dif- 
ferent sub-interval of (e,r<,(e)). These functions match 
along the overlaps between the sub-intervals. In terms of 
the variable r, they read 
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zLir) , „ „ + , , 9 +••• ; inner region : e < r < \/e (3.29) 

el 1 e'' 3 

^™<i(^) = ^ + ~ + 2^'^ + 2 ^ ^ 2^^^^ ' middle region : e < r < 1 (3.30) 

r 1 

• • ; outer region : ^/e <^ r << 1 - \/e (3.31) 



z' 

equ 



^/j— 72 ^2)3/2 

1 + ie 

^ + • • • ; equatorial region : 1 — r ^ 1. (3.32) 
V2(l -r~e) 



2. Longitudinal stretching 

In order to get the actual profile z{r), one has to 
integrate equations (|3.29I3.32|I . Since the four regions 
overlap each other along finite intervals, artificial lim- 
its of integration (within the overlapped intervals) must 
be introduced. Thus, we integrate (|3.29l) from e to e" 
(1/2 < a < 1), (j^^Ujl from e" to e'' (0 < (3 < 1/2), 
(|S3T|) from e'' to 1 - e'' (0 < 7 < 1/2), and from 
1 — e''' to 1 — e. As will be verified below, the values of 
a, /?, 7 do not appear in the final expressions to leading 
order in e. 

Performing the integrals of equations H3. 2913. 321) be- 
tween the above limits, we find the leading contributions 
to be: 

/i y z[„dr w (a- 1) eloge-t- • • • 
1 

h = / zL,dr«l- V2e^/2- -6^^-/3 e log £+••• 

Jel^ 2 

I4 = I z^,„dr « V2e''/2 _^ • • • (3.33) 



l-£7 



In particular, the z-coordinates of the centers of the 
spheres (in practice, the 'equatorial' values of z at which 
dz/dr — 0) are shifted from ±1 (for e = 0) to ztzg ~ 
±z(r,(e)): 

=/i +/2 +/3 + /4 ~ l + e|loge|. (3.34) 



B. Dynamics 



In a recent work ijDavidovitch et oil . l2002(l we proved 
that all CMC surfaces, except the plane, are unstable 
under the dynamics H2.7|l . This holds in particular for the 
fixed point surfaces described by equations (|3.29I3.32|I . 
Consider a surface x close to a fixed point surface x*: 



x(x*) =x* +^(x*)n(x*). 



(3.35) 



where n is a unit vector normal to the unperturbed sur- 
face, and S is the magnitude of the perturbation. Linear 



stability analysis of the dynamics H2.7|l near x* results in 
the formal equation for the perturbation 6: 



(3.36) 



where £ is a linear differential, generally non-Hermitian, 
operator, whose coefficients depend on the geometry of 
the surface x*. In this section we will compute, to 
leading order in e, the largest eigenvalue Xmax{^), for 
X* a CSCMC whose profile z{r) is described by equa- 
tions (|3.29I3.32|I . Xmax determines the rate of growth 
of a generic perturbation S with a non-zero component 
along the maximal eigenvalue of equation (|3.36(l , whether 
it leads to sintering of the spheres or to dissolution of the 
neck. 

Derivation of the detailed form of C requires some tools 
of differential geometry, and will not be presented here. 
For o ur purposes it is enough to use the following re- 
sult (Davidovitch et "oZI . i2002) : The eigenvalues of C are 
identical to the eigenvalues of a Hermitian operator TL: 



(3.37) 



where is a surface gradient, and the 'potential' V^(x*) 
is: 

y(x*) = 2Hl-K~]^[VlB + {WsBf] 



B(x*) 
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l0g(i/2-if(x*)), 



(3.38) 



where i?,, K = KaKi are respectively the mean and Gaus- 
sian curvatures of the unperturbed CMC surface. In par- 
ticular, Xmax{^) = XmaxiTi-)- For a CSCMC with profile 
z{r) we use the coordinate system (f, z, (p) to write: 



K{r) 



1 



1 + z' 



(1 + Z'2)2 



r(l + z'^) 



/2^2 ■ 



(3.39) 



We note that for the CSCMC described by equa- 
tions H3. 2913. 32(1 . the Gaussian curvature K{r) and its 
derivatives diverge toward the neck (r — > e), as e — > 
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0. This can be easily understood by recalling that 
K — Kai^i- The largest value of the azimuthal curva- 
ture Kq = l/r is clearly achieved at the neck, where 
Ka,neck = l/^- Since for this surface k; = 2 — Kq « — 1/e, 
we see that Kneck ~ — as e ^ 0. As one moves 
away from the neck to the spherically shaped periphery, 
Ka ^ 1, the Gaussian curvature K ^ 1, and therefore 
\V\ decreases. Thus, the maximal eigenvalue XmaxO^) is 
determined by the region near the neck, and we can com- 
pute it by considering the operator H. in the inner region 
e < r <C -s/e only. According to equations H3.2()I8.39|1 , the 
change of variables: 



(3.40) 



enables us to write = 9^ in the inner region. Some 
tedious algebra (see appendix) yields the potential form: 



V{x) 



3 e^-l-ix^ 
4(e2 -f- j;2)2^ 



(3.41) 



in the inner region. With the Id potential (|3.41ll . the 
eigenvalue equ ation for Ti is equivalent to a s pheroidal 
wave equation ijAbramowitz and Stegunl . Il97(tl . Numer- 
ical calculation of Xmaxi'H) in the limit e — s- yielded 
A™a.(W) = 0.282e-2 + o(£-i). 

We conclude that the profile given by equations H3.29t 
13.32(1 is an unstable fixed point of the dynamics H2.13|l . 
Typical perturbations will increase at a rate proportional 
to e^^, and will lead either to sintering of the two spheres, 
or to dissolution of the neck. 

This has some immediate consequences. Consider, for 
example, two spheres compressed against one another so 
that the distance between their centers is d < 2R. A sin- 
ter formed around the contact zone between such spheres 
will generally grow, since the contact geometry should 
typically be on the growth side of the critical neck geom- 
etry. On the other hand, if two spheres become connected 
by a neck while their centers are at a distance d > 2R, 
the dynamics of the neck may cause it to either grow 
or evaporate, the boundary between these two behaviors 
being approximated by d « 2i?(l + e\ loge|). 



IV. ASYMPTOTIC DYNAMICS 

In this section, we will study the asymptotic behav- 
ior of two sintering spheres when the surface profile is 
far from any fixed point, but the neck size is still small 
compared to the radius of the spheres. In this regime, 
the non-linear dynamics l|2.13|l cannot be approximated 
by its linearized form (|3.36() . In figure 3 we show the 
numerically computed evolution of an asymptotic neck 
profile, which appears to evolve according to a similarity 
law. 




FIG. 3 Evolution of cylindrically symmetric sintering profile 
under equation 12.1311 . The initial profile is parabolic. 



For volu me preserving d ynamics similar to l(2.13|l . 
Kuczinsky l|Kuczinskvl . Il949l) suggested that the asymp- 
totic neck profile evolves like a spherical cup, 



r(z, t) = rait) + w{t) - ^/w^{t) - z^; 0<z<w{t), 

(4.1) 

where ro(t) is the neck thickness and w(t) its width. This 
neck profile has to be patched to the static peripheral 
spherical shape: 



"{z,t) = \j2z - z2 



(4.2) 



Patching the z-coordinates of the two regions at z ~ 
w{t) implies: 

{w{t) + rQ{t)f « 2w{t) w(t) « ra{tf/2 < ro(i), 

(4.3) 

where we used ro(i) <C 1. To obtain the time dependence 
of rQ(t),w{t) we recall that the RHS of equation ((2.13|l is 
proportional to \{K.a + ni) ~ 1, where Ka,neck ~ l/?'o(i), 
and the spherical cup profile (|4.1() implies ki « —l/w{t). 
According to equation ((4.3(1 w{t) <C ro(t) <C 1, and there- 
fore the RHS of ((2.13(1 is dominated by w{t)~ , whereas 
the LHS is simply dro/dt. Using equation 1(4.3(1 we get 



drp 
dt 



1/3 



(4.4) 



in agreement with experiments. A similar derivation for 
sintering dominated by surface diffusion yields rg ^ i^^^, 
and other gro wth rates can be derived for other kinds of 
kinetic routes f German, '1986*; 'H errind. ll95Ct iKuczinskvl 
^M9; Swinkcls and Ashbv, 1981). 

It should be stressed here that Kuczinsky's profile is 
simply a guess and is not in any sense an asymptotic solu- 
tion of equation ((2.13() or of its volume preserving version. 
Moreover, notice that the t^^^ growth rate cannot be de- 
rived for a general neck profile, but rather depends on 
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the assumed shape. As an example, consider a paraboHc 
neck profile: 



r{z,t) = r{t)+a{t)z^. 



(4.5) 

-> of 



Patching profiles and derivatives with the limit z 
(|^ gives ro - t^/^. 

Kuczinsky's assumption of semi-circular neck profile 
(|4.1|l was questioned by several authors, and other pos- 
sible asymptotic neck profiles were suggested for sin- 
tering dominated by evaporation-precipitation and by 
other kinetic route s . In particular, sev e ral re searchers 
ijAmar et all Il989t iGerman and Muniit Il975j) . conjec- 
tured that in a variety of sintering processes, the neck 
profile evolves as a CSCMC surface, whose spatially con- 
stant mean curvature evolves in time. Notice the simi- 
larity between this assumption and Kuczinsky's one, ac- 
cording to which the longitudinal curvature ki is spa- 
tially constant at the neck, whereas Ka may slightly 
vary. Other researchers, who focused their analysis 
on sintering dominated by surface diffusion, performed 
careful numerical simulations and advanced analytic ar- 
guments, from which they concluded that the asymp- 
totic neck shape evolves very differently f rom the shape 
suggested by Kuczinsky, equation 114. ID. llEgg ers''. 1998; 
iGerman and Lathroi)lll978HNichols and Mullina. il965i|) . 
Moreover, it was shown by these researchers that in that 
case the asymptotic growth rate of the neck is not t^^"^ , 
as was suggested by Kuczinsky-like arguments for sur- 
face diffusion dynamics, but rather becomes closer to 
^1/6 Qj. gven t^/^. Other studies, which focused on ki- 
netics dominated by surface and grain boundary diffu- 
sion showed that the neck profile and its growth rate in 
this case are also much more complicated than t he sim- 
ple models suggested by Kuczin sky (Bross and Exxneil 
ll975lSwiTike]s a,Tid AshbvL 1198(1. 

Here we show that the asymptotic neck profile for the 
sintering process determined by equation H2.13|l evolves 
to a similarity shape, which is indeed different from 
Kuczinsky's or other suggested shapes. We show how- 
ever, that Kuczinsky's t^^^ growth rate is correct in this 



Let us write the sinter profile r{z,t) as 
r{z,t) ^roit) + b{z,t), 



(4.6) 



with 6(0) = 0, so that ro{t) is simply the radius at the 
center of the sinter. We suppose that throughout the 
sintering, non-spherical, region, we can assume 



b{z,t) < roit), 



(4.7) 



an assumption we must verify for self-consistency at the 
conclusion of the calculation. Then we have 



b" 1 

^ w 2ro H = a(t), 

1 + b'^ ro ^ 



(4.8) 



where we have replaced r by Tq and 2yl + r^ « 0, their 
values at z = 0, while retaining the longitudinal curva- 



ture term (which also dominates the curvature in Kuczin- 
sky's solution). The solution of equation 14.8|l is: 



6(z,t) 



a{t) 



log cos[a{t)z\, 



(4.9) 



where we assumed the boundary condition: fo'(O) ~ 
(from z ^ —z symmetry). 

Following a similar reasoning to the one presented in 
the beginning of this section, we notice that patching 
derivatives of the concave neck profile given by equa- 
tions (|4.6I4.9() , and the convex spherical periphery, equa- 
tion (|4.2I) is possible only for (r, z) — *■ (r*, z*) where both 
derivatives become much larger than one: 



dr 



from which we get: 



2a{t) 



+ 0(a-3/2). 



(4.10) 



(4.11) 



Matching to the peripheral spherical shape equation H4.2|) 
implies: 



From equations (|4. 814. 1114. 12|l we get: 



1 TT TT 

ro 



(4.12) 



(4.13) 



whose asymptotic solution [at ro{t = 0) ^ ro(t) <C 1] 
gives the Kuczinsky growth rate (|4.4(l , 



ro{t) « ( — 



1/3 



(4.14) 



To check the consistency of this solution we must ver- 
ify the validity of the approximation |6(z, i)/ro(i)| ^ 1 
for z < z* . This seems clear, since although b' is large at 
the patching point for small z* , b is still small (6 <C tq) 
because of equations (|4.9I4.11|I . However, while inside 

the neck, near z = 0, the term 2\/ 1 + z'^ (excluded from 
equation H4.8|) ')is negligible, we should insure that it does 
not become large as z — > z*. Therefore this term should 
be included in the consistency check, and we seek a so- 
lution to the equation: 



b" 



2Vl + b'^ = 2ro + — = 



1 + 5,2 . ^v. , . -Zro + -^ait), (4.15) 

under the same BC as above: 6(0) = 6'(0) = 0. With the 
change of variables: tanu = b', equation (|4.15|) recasts 
into the form: 



u' + 2yr+~tai?u = a{t), 



(4.16) 
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from which we get 



References 



b{z*) = [ dzb' = f 
Jo Jo 



tan"^[6'(2*)] 



du- 



tan u 



'.{t) - 2%/TTtai? 



2a{t) 



log(l + [5'(z*)]^ 



(4.17) 



This result is identical for small z* to equation 1)4. 9|l . from 
which we verify the consistency of equation ()4.1(J|) . 



V. CONCLUSIONS 

In this paper we studied a sintering model that we be- 
lieve to be relevant for understanding the evolution of 
morphologies in sedimentary rocks and other geophys- 
ical systems, under reaction-limited kinetics and open 
environmental conditions. The assumption of reaction 
as opposed to transport limited kinetics is quite de- 
fensible in this context. For example, for siliciclastic 
rocks (quartz in water) we take as representative val- 
ues: L 



par 



^pore 



10-'*m, D « 10-'-*m7s, and 
Kf 2 lO^^^m/s. From equation H2.6|l we obtain the 
typical time scales Tdiss ~ lO^^s » Tdiff ~ 10s. However, 
much of the diagenesis of sedimentary rocks is driven 
by pr essure so l ution , rather than b y surface energy ef- 
fects (Ide Boeii 119771 iRutted . \l97(ti . Surface energy ef- 
fects might dominate in systems in which the solid phase 
has not experienced significant pressure changes over the 
lifetime of the grains-a situation that might prevail in 
magmatic plumes, in which solidified material coexists 
with liquid magm a still satur ated with the components 
of this material (Visseilfl999fl . 

In glaciology, the most likely complication to our as- 
sumption of dissolution-precipitation kinetics driven by 
surface energies is likely to be surface melting, especially 
near the melting temperature of the ice. Transport via 
surface diffusion or in melted surface layers is also likely 
to be significant for sintered metallic powders, especially 
if the vapor pressure of the metals is low. 

Even though the sintering physics of real technological 
and geological materials is likely rarely to be in a regime 
where the rather simplified model of this work includes 
all of the relevant physics, we nevertheless believe that 
our methods, based on careful asymptotic analysis of the 
stationary states of the kinetics, combined with similar- 
ity analysis of late-stage sintering, not only illuminate the 
physical limit we have chosen, but also should deal suc- 
cessfully with other types of tranport mechanisms. This 
is likely to be the subject of subsequent works. 
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APPENDIX 

In this appendix we prove that the potential V, as 
defined in equation (|3.38() has the form H3.41|l in the inner 
region Zi„(r). For the CSCMC surfa ce t\:2i)l\.:]'2\i H ^ 1, 
WsH = 0, and from equation 1)3.38(1 we get: 



V = 2-K- 



1 V^A' 9 (VsK)' 



mi-K 128(l-is:)2 



(5.1) 



From equation H3.39|l we get: 
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1 



r(l + z' 
1 



/2-)5/2 



z"' + z'z'" ■ 



1 + Z' 



/2 



r(l 



.'2^3 



[ 3zV" + zV"' + 2 



z'z'" , z"^ 



z'z"^ z'^z"z"' 
13 ^ - 13- 



1 + z' 



.'2 



28- 



9- 



/2 //2 

Z Z 



(1 + Z'^)2 ■ >(1+Z'2) 



(5.2) 



Consider now the profile Zm(r) in the inner region e < 
r <C v^, as given in equation l|3.29|l . Since in this region 
r^/(r2 — e^) ^ .^^g gj^j ^j^at the leading contributions 
to the profile derivatives are: 



Substituting the last expressions for the profile deriva- 
tives in equation (|5.2I) . and using the change of vari- 
ables (|3.40() . we obtain equation (|3.41() . 



(r 



2 _ g2^1/2 ' 



er 



£2)3/2' 



(r2 - e2 
6er^ — 9re^ 



(r2 



2)5/2' 



(r2 



2)7/2 -(S-S) 



